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The flow equation method (Wegner 1994) is used as con- 
tinuous unitary transformation to construct perturbatively ef- 
fective Hamiltonians. The method is illustrated in detail for 
dimerized and frustrated antiferromagnetic S = 1/2 chains. 
The effective Hamiltonians conserve the number of elementary 
excitations which are 5 = 1 magnons for the dimerized chains. 
The sectors of different number of excitations are clearly sep- 
arated. Easy-to-use results for the gap, the dispersion and 
the ground state energies of the chains are provided. 

PACS numbers: 75.10.Jm, 02.30. Mv, 03.65.-w 



I. INTRODUCTION 

Perturbation theory is one of the most important and 
most versatile tools for problems which are not exactly 
solvable. Various methods depending on the problem 
under study have been invented and used. Due to the 
enormous increase in computer capacity it is a very in- 
teresting task to use algebraic programmes to perform 
perturbative calculations. 

The aim of the present work is to propose a general 
perturbation scheme which splits naturally into two sub- 
sequent steps. Both these steps can be implemented in 
a direct manner on the computer. The first step is not 
model specific. It relies only on two prerequisites. 

(i) The unperturbed Hamiltonian Hq must have an 
equidistant spectrum bounded from below. With- 
out loss of generality we may assume that Ei = i 
for i e {0, 1, 2,3,.. .}. We say that i denotes the 
number of energy quanta in the system. By Ui the 
corresponding subspaces are denoted. 

(ii) The perturbing Hamiltonian iJs links subspaces Ui 
and Uj only if \i — j\ is bounded from above, i.e. 
there is a number iV > such that Ha, can be 
written as iJs = Sn=-A' ^n where r„ increments 
(or decrements, if n < 0) the number of energy 
quanta by n 



[ifcTn] = nTn 



(1) 



Thus the full problem reads 



N 



H = Ho + X J2 ^" 



(2) 



~N 



where A is the perturbation parameter supposed to be 
small A < 1. In this work we will restrict to A^ = 2. 
The first step consists in finding a systematic mapping 
of the problem in Eq. (0) to an effective one given by a 
Hamiltonian Hcs which conserves the number of energy 
quanta. 

The second step is the model specific one. It consists 
in the actual calculation of H^s for a given number of 
energy quanta. 

To illustrate the abstract ideas formulated above we 
will use the frustrated and dimerized 5=1/2 chain given 
by (j counts the sites) 

L 

H = JoY,[{l + i-iyS)S,S,+, + aoS,-iS,+i] , (3) 

where L is the number of sites. A situation consistent 
with Eq. (y) is found for strong dimerization. Hence 
we rewrite Hamiltonian (H) as (subscript i counts the 
dimers) 



H — J 2_^ [S2iS2i+l + AS2iS2i-l + 

A a(S2iS2i-2 + S2j-iS2j+i)] (4) 



i=0 



with 



j^Mi + s) 

X = il-S)/il + S) 
a = ao/(l -S) . 



(5a) 
(5b) 
(5c) 



The unperturbed part Hq (up to a trivial constant 3L/S) 
and the perturbing part Hg are then 



Ho = J2 [S2«S2»+i + 3/4] 



(6a) 



Hs — 2_^ [S2iS2i-l 



Q;(S2iS2i-2 + S2i-iS2i-i-i)] , (6b) 
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where we measure implicitly all energies in units of J. 
The ground state of Hq is the product of singlets on the 
dimers, i.e. the bonds (2i, 2i-\- 1). The energy quanta are 
here the excited dimers, namely the local triplets. The 



number of triplets classifies the degenerate energy eigen 
spaces of the unperturbed problem. 

Besides the purpose to serve as an example for per- 
turbation by flow equations the frustrated and dimer- 
ized chain is of considerable physical interest itself. Ideal 
spin-Peierls systems are one-dimensional spin systems 
which are coupled to the lattice. At low enough tem- 
peratures they dimerize since this dimerization leads to 
a gain in magnetic energy oc 5^'^ which overcompen- 
sates the loss in elastic energy oc 5"^, see e.g. ||l[ and 
refs. therein. So spin-Peierls systems provide dimerized 
spin chains in a natural way. The first inorganic spin- 
Peierls substance CuGeOa in particular provides the ex- 
ample of a frustrated and dimerized spin chain since 
there is much evidence that a certain amount of frustra- 
tion is present in this substance [0-0|. Other examples 
are strongly anisotropic substances where the dimeriza- 
tion is built-in in the chemical structure. Examples are 
Cu2(C5Hi2N2)2Cl4 i and (VO)2P207 §. Of course, 
the real substance mostly display also some additional 
two- or three-dimensional coupling. But the approach 
we present here is suited to tackle even these systems, 
see e.g. Ref. |j. 

One might argue that exact diagonalization or quan- 
tum Monte Carlo approaches are better suited to calcu- 
late dispersions u}{k) or similar quantities in d = 1. These 
methods, however, yield only the result for the chosen pa- 
rameter set. The perturbative results, however, will be 
obtained as polynomials in the weak bond coupling A and 
the frustration a. Thus, once computed, anyone can use 
the perturbative results easily to fit measured or other- 
wise obtained data. Thereby an extremely fast method 
for the determination of coupling constants is provided. 
Of course, the perturbative approach can be applied only 
for A < 1 where the equal sign represents the worst case. 
For A > 1 the perturbative approach breaks down. 

The work is organized as follows. In the next section we 
extend the approach of Stein [g[ who did a calculation for 
iV = 1 up to fifth order to A^ = 2 and up to tenth order. 
The work of Stein improved earlier calculations^] which 
generated more intermediate terms (see Ref. ||§] for dis- 
cussion). The flow equation transformation which is used 
by Stein and by us was introduced by Wegner five years 
ago |lQ]. In Sect. Ill we illustrate our method by apply- 
ing it to a one-dimensional Heisenberg antiferromagnetic 
5 = 1/2 chain. The effective one-triplet Hamiltonian is 
computed. The ground state energy, gaps and dispersion 
relations are discussed in Sect. IV. Summary and outlook 
conclude the main part of our work. 



II. PERTURBATION BY FLOW EQUATIONS 

All what is presented in this section is based only on 
the fact that the initial problem has the form (0) fulfilling 
the requirements (i) and (ii) with N = 2. 

The general idea behind the fiow equation approach 



introduced by Wegner |10| is to perform a continuous 
unitary transformation which makes the problem more 
easily tractable. Mostly, one tries to make the problem 
"more diagonal". In our case we will achieve a block 
diagonal form. A broad field of application is to iden- 
tify certain quasi-particles for which an effective Hamil- 
tonian can be found. Here, flow equations can be used 
to implement a renormalization of a given problem on 
the Hamiltonian level, not only on the level of certain 
observables or couplings |10|-[l^. Analogous ideas were 
suggested parallely by Glazek and Wilson in the form of 
similarity transformations p3| . 

In the present work we do not focus on the renormaliza- 
tion properties of the flow equation approach. Following 
Stein 1^, we use them to implement in a systematic way 
a continuous unitary transformation which maps the per- 
turbed system onto the unperturbed one which is easy to 
understand. 



A. General Formalism 

According to the original idea a running variable i is 
introduced which parameterizes the continuously evolv- 
ing Hamiltonian H{£). The starting operator is the bare 
Hamiltonian; the operator at inflnity is the desired effec- 
tive Hamiltonian 



H{0) 

H{oo) 



HcS 



XiT^2+T^i+To+Ti+T2 



(7a) 
(7b) 



The unitary evolution is engendered by its antihermitcan 
infinitesimal generator ri{£) 



dH{£) 
di 



^W),H{£)] 



(8) 



Applying naively Wegner's choice for the generator 
ri{i) = [Ho,H{£)] the resulting differential equations 
quickly become very messy since the band block diag- 
onal structure of the original problem is lost. By "band 
block diagonal" we mean the fact that A^ has a finite 
value which does not change in the course of the flow 

We will choose a slightly different infinitesimal gener- 
ator which allows to keep the band block diagonal struc- 
ture of the original problem, i.e. the parameter N stays 
2 for all values of £. The general Hamiltonian H{£) can 
be written as 



H{£) = Ho + Xe{£) 



(9) 



and the operator Q{£) links only subspaces Ui and Uj 
with \i-j\ < 2. 

The most general form of O(^) is 



e(^) = ^A^-i J2 ne;m)T{m) 

|m|— fc 



(10) 



fc=i 



where F{t, m) are real- valued functions for which we will 
derive nonlinear, but recursive, differential equations be- 
low. The other symbols are 



771 = {nil, 7712,7713, 

771, e {0,±1,±2} 
|to| = k 

1 \7Tt) -^nii-^ni2-^rii3 • 



, TOfe) with 



■ T„. 



M{m) = Y.' 



(11a) 
(lib) 
(lie) 
(lid) 

(lie) 



i=l 



The vector m together with the product T{rn) encode 
all possible products of the incrementing (decrementing) 
operators T„ as defined in Eqs. (1^,0) ■ The infinitesimal 
generator of our choice reads 



77(£)=^A''' Y^ sgn{M{m))F{£;m)T{r_ 



(12) 



fe=i 



This choice is very similar to the infinitesimal genera- 
tor Mielke proposes for band matrices ||l^. An adapted 
version of his proof that such an TTleads to (block) diag- 
onality can be found in appendix H. For the purposes of 
the present perturbative approach the general proof can 
be replaced by the observation that the transformation 
can be performed successfully to all finite orders. This 
will be shown below. 

A short computation shows that substituting 
sgn(M(m)) in Eq. (|2|) by M{rn) would correspond to 
the first choice ■7i{£) ~ [Ho,H{£)]. This relies on 



[Ho,T{m)] ^ M{m)T{rn) 



(13) 



which is a straightforward generalization of Eq. (|l|). 

Insertion of the ansatze nu) and (n3) into Eq. (H) leads 
to 

A^ = A[ry(£),e(^)]- 

oo 

^A*-' Y^ sgn{M{m))F{£;m)[Ho,T{m)] . (14) 

k=l \m\=k 

Comparison of the coefficients for each term T{rn) yields 
then a differential equation for the functions F(£; 7n) 



Pii;in.) 



-\M{m)\F{£;m)- 



(15) 



J2 [sgn{M{mJ) - sgn(M(m2))] F{£;rn,)Fi£;m2). 

The summation condition {tu^, TOqI — jn means that one 
sums over all possible nontrivial breakups of m 

ni l — ("^l) ^-Ild TO2 — ('7^2, ■ ■ ■ , TTlk) 

nil ~ ('^ii '712) and m^ — {m^, . . . , m^) 



nil ~ ('^ii '^2, 'TT-s) and TO2 — {rrii. 



,mk) 



nil = ('7^l, • • • , rnk-2) and ni2 = (mfc_i, m^) 

mi = (mi, . . . , mfe_i) and ni2 = (mfe) . (16) 

This summation notation will also be used in the follow- 
ing. The starting conditions follow from ( |7a| ) 



F{{);m) — 1 for |m| — 1 
i^(0;m) = for H>1 



(17a) 
(17b) 



From Eqs. (HS 17) we can deduce a number of relations 
by induction. First, we see that the functions F{£;ni) are 
always real. Furthermore, they obey the two symmetry 
relations 



F{£;-m) ^F{£;m) 
F{£;-rn) = (-l)l2il+ii^(£;TO) 
where we use the notation 



m 



{mk,nik-i, ■ ■ .,TO2,mi) 



(18a) 
(18b) 

(19) 



The square bracket in Eq. (y_5|) ensures that the sum van- 
ishes if I M(m) I > 2 



F{£;m)^Q for |Af(TO)| > 2 



(20) 



For instance, a term generating three energy quanta 
M{inri) — 3 could only be induced from terms with 
M{rrii) = 2 and M{ni2) = 1 or vice- versa. But such 
combinations are suppressed by the square bracket in 
Eq. ([1^) . This observation is at the basis of the preserva- 
tion of the band block structure [ 14 1 . If we had chosen the 
infinitesimal generator ri{£) in (12) without the signum 
as it would correspond to Wegner's original sug gestion 
ri{£) = [Ho,H(£)] the square bracket in Eq. ( [l^ ) would 
read [M{nii) — M{rri2)] and hence the band structure of 
the couplings would be destroyed for £ > 0. 

For the solution of Eq. (|l^ we observe that the first 
term on the right hand side just generates an exponential 
prefactor 



F{£;m) = exp{-\M{rri)\£)f{£;m) 



(21) 



The rest of the equation dl3) is recursive and can thus be 
directly found by integration beginning from the starting 
conditions 



d 

d£ 



f{£;m) = V e(l*^(2i)|-|A^(™i)l-|A^(™2)l)' . 
[sgn{Mimi))-sgniMini2))]f{£;mi)fi£;m2) . (22) 



Note that \M{ni)\ - \M{nii)\ - \M{m2)\ < holds so 
that no exponential growth occurs in the f{£;7n). Let us 
focus on the functional form of the functions f{£; m). We 



state by induction that sums of terms te "^^^ occur for 
non-negative integers i and /x. More precisely, we obtain 



r(m) 



(23) 



where the degree of the polynomials P^ (i; m) is always 
equal or less than \rn\ and the upper limit of the sum 
r(TO) obeys 



r(m) = 2 



if \rn\ \ 

- -|M(m)|+^|m,| 



(24) 



In principle it is also possible to write down explicit re- 
cursion relations for the polynomials, see for instance Ref. 
pi for iV = 1. But they are of little clarity. If the actual 
calculation is done by symbolic calculation it is sufficient 
to retain that according to Eq. dlS) or Eq. ( p2[ ) expres- 
sions of the type ( p3| ) have to be multiplied, added and 
integrated. This is a straightforward task and can be 
implemented in symbolic programmes. 

The quantities we are finally interested in are the coef- 
ficients of iJcff = H{i = oo). From Eq. ( pi] ) we know that 
only terms with M{rn) = will not vanish for ^ ^ cxd. 
This is exactly what we intended to achieve since we want 
ffoff to commute with Hq (cf. Eq. (13)) so that the num- 
ber of energy quanta (triplets in our example) becomes a 
conserved quantity. Hence, we can write the final result 



H,g = Ho + J2^'' Yl C{m)T{m) (25a) 

fc=l |m|=fc,J\/(m)=0 

C(m) = F(oo;m) . (25b) 

Further details on the computation of the coefficients will 
be given in the next subsection. Results for the C(m) are 
presented in appendix pi 



B. Computer Aided Evaluation 

We implemented the coefficient computation in C-\ — h 
because of its high performance and its class concept 
which we used to encode the basic data elements. As 
an example let us consider a generic fourth order coeffi- 
cient (i.e. k — A) f{l;rn) : 



1 



/(Z;(l,l,l,-2)) = ^; + ^e-2'--e-4' 



2 2 
This can be stored as a list of basic data elements like 



(26) 






(27) 



each containing four separate integers p, q, i, ^. In fact, p 
and q can become very large. So they have to be stored 



in a multiprecision data type like long long int on some 
Unix systems. Still, all computations can be done very 
fast in integers and the results are rigourously exact. 

Eq. ( p2|) is essential in the computation of the f{l;m). 
The basic idea is to build two loops. The outer loop 
controls the order starting at k — 2 since the initial con- 
ditions are the values of the f{l',m) in first order. The 
inner loop generates all possible m in the current order 
k. A single rrii in m ~ (rrti , . . . , mi, ... , ruk) can take five 
different values, see Eq. (lib). Thus we can introduce 
the loop variable n G {0, 1, . . . S'^ — 1} 



E«»5^ 



(28) 



j=0 



The coefficients a^ G {0,1,2,3,4} are mapped uniquely 
onto the set {—2, —1, 0, 1, 2}. Finally we retain those to, 
with |M(to)| < 2. In this way, Eq. ( p^ ) can be calculated 
for each k and m. 



The calculation of a single f{l;m) via Eq. (22) can be 
spht in four steps, (i) One has to encapsulate Eq. ( [2^ ) 
in yet another loop controlling all possible breakups of 
TO (cf. Eq. dig)), (ii) The functions f{l;mi) and f(l\m2) 
known from calculations in lower orders have to be mul- 
tiplied for each breakup of to. (iii) One has to sum over 
all breakups, (iv) Finally, the result from steps (i) to (iii) 
has to be integrated. 

Since p and q in the basic data elements (E7h can be- 
come very large during addition and multiplication, both 
operations employ Euklid's algorithm to generate maxi- 
mally canceled fractions p/ q. To minimize memory us- 
age these operations contain simplification subroutines 
based on the quick-sort algorithm. These subroutines 
sort according to increasing powers of I and e^^'. Simul- 
taneously, addends are identified and added if they are 
of equal type. The resulting f{l;m) consists of linearly 
independent addends only. Due to the quick-sort algo- 
rithm the computation time as function of the number of 
addends n is only of order nln(n). 

The final integration can be done easily. The functions 
f{l;rri) break down to basic data elements ( p7| ) so that 
(a>0) 



dl'l" = , 

» + l 



i+l 



(29a) 



/' dr/'^e-"'' = - . J 1 - e-"' y -^— \ (29b) 

achieves the integration. 

To calculate the C(to) one has to perforin the I — » oo 
limit on those f{l;rn) for which \M{rn)\ — 0. Note 
that for |M(?7i)| = one always has a > since 
a = |Af(TO)| - \M{rn^)\ - |M(to2)| = and \M(m)\ = 
implies the vanishing of |M(to-|)| and of \M{rri2)\- Hence 
the right hand side of Eq. (|2^) vanishes due to the 
square bracket containing the sign functions. So no such 



f{l;rn) is ge nerated. Thus the C{rn) can be calculated 
by Eq. ( 29b ) omitting the term proportional to e~"' on 
the right hand side. 

The symmetry relations (|l^) can be used as a check of 
the results. The C{rn) are saved in a file together with 
the corresponding m for later usage. They are given up 
to order fc = 6 in appendix H. We intend to provide them 
up to order fc = 10 in electronic form on our homepages 
on appearance of this article. 

Unfortunately, Eq. (22) implies also a natural limita- 
tion of the computation. Because of its recursive nature, 
the f{l;rn) of all preceding orders have to be stored. 
They are needed to derive the f{l;m) in the current 
order. This leads to an exponential memory increase. 
To calculate all C{rn) to order fc = 10 inclusively, we 
used about 30,000,000 basic data elements ( |27| ) occupy- 
ing about 1GB RAM. Because of the extensive memory 
use we employed a SUN Ultra Enterprise 10000 which the 
Regional Computing Center of the University of Cologne 
kindly placed at our disposal. The calculation took about 
12h. 



III. APPLICATION: 
DIMERIZED AND FRUSTRATED S = 1/2 CHAIN 

In this section we demonstrate how the knowledge of 
the C{rn) in the effective Hamiltonian i/off ( pq ) permits 
to perform specific calculations. The first step is to eval- 
uate the operators r„ for the model under study, here 
dimerized spin chains. Then we calculate the ground 
state energy, the energy gap and the one-magnon disper- 
sion of the one-dimensional Heisenberg antiferromagnetic 



SoiSo 



T_ 



-T_i+To+Ti+7-2 



(33) 



S 



chain. This is again done by implementing the 



requiring that the subscript indicates the net change of 
the number of triplets. In other words matrix elements 
connecting a ket of two singlets with a bra of two triplets 
belong to T2 and those connecting a ket of one singlet 
and one triplet with a bra of two triplets belong to Ti 
and so on. In this way one finds all the T„ and their 
sum is S2iS2i-i. Table m summarizes the results. For 
later convenience we split 70 = 7^+7^. All other ma- 
trix elements can be constructed by using the relation 
TJj — T_„. To incorporate the effect of frustration it is 
sufficient to note that a triplet is invariant under spin 
exchange whereas a singlet acquires a factor -1. Let 



a{S2iS2i^ 


-2 + S2.- 


1^2i+l) 


= y^^n- (34) 
n=-2 


47S 


\t''^\s) 


— >■ 




-|s,i"'±^) 


47S 




■ > 

■ > 

■ > 

■ > 




|tO,i»>-|f±\iTi) 

\t\t-^) + \t-\t^) 


40-1 


\s,t'},\t',s) 

\s,t'>),\t°,s) 

\s,t~''),\t'\s) 


>■ 

> 

y 




\t\t-^)-\t-\t') 
\t'\t-')-\t-\t") 


4T2 


\s,s) 


■ >■ 


\t' 


i-^>-|t",i") + |t-\i^) 



equations on computer. 



TABLE I. Action of the operators 7i as defined by Eq. (|3^ 



A. General Equations 

The explicit form of the operators r„ has to be deter- 
mined so that 



Hs=T^2+T^i+To + Ti+T2 



(30) 



(cf. Eqs. (p 6b). Let us consider one addend of H^ as 
starting point 



S2jS2i-l + a (S2iS2i-2 + S2i-lS2i+l) 



(31) 



Obviously, only neighbouring dimers are affected. For 
simplicity we first calculate the matrix elements 



(a^i-i, Xi\S2iS2i^i\xi-i,Xi) 



(32) 



where Xi-i,Xi G {s,t^,t°,t~^} are singlets or one of the 
triplets on the adjacent dimers z — 1 and i. The super- 
script n e {0, ±1} in t" stands for the S^ component. 
For some fixed value of i we write 



By spin exchange the T'j operators defined in Eq. (|3^ 
are reduced to the T,; as given in Table | 



T^2 ^ ^2a • 1^-2 



T- = 2a . T« . 



Finally, Eq. ( pO|) implies 



T±2 =J2il-2a)- T±2 

i=0 
■§■-1 

i=0 



To = E (1 - 2a) • Tg + (1 + 2a) • Tg 



i=0 



(35a) 
(35b) 
(35c) 
(35d) 



(36a) 
(36b) 
(36c) 



The subsequent subsection shows how the r„ operators 
can be implemented. 



For the antiferromagnetic S — ^ Heisenberg chain 
given by Haniiltonian (|6|) we can then calculate the 
ground state energy by 



Eo = (0|i/eff|0) 



(37) 



where we used the shorthand |0) :— |J|j Si) for the prod- 
uct state of singlets on all dimers which can be viewed as 
triplet vacuum. 

The actual calculations are done on finite clusters. 
From the linked cluster theorem we know that the finite 
order contribution of a short-ranged perturbation does 
not depend on the cluster size for sufficiently large clus- 
ters. In our one dimensional example only neighbouring 
dimers are linked. So it is sufficient to consider II dimers 
to avoid a wrap-around in order 10, i.e. for II dimers 
or more we are sure to find the thermodynamic contri- 
bution. Moreover, we can check the size-independence 
explicitly. 

To calculate the one magnon dispersion we have to 
consider the subspace with exactly one single singlet be- 
ing excited to a triplet. Using |j) = |s, s, . . . , i, . . . , s) 
with one triplet t (S^ component does not matter) on 
dimer j we compute the action of effective Hamiltonian 
(Eq) on \j). Since the number of triplets is conserved by 
construction the initial triplet can only be shifted. 



Hcslj) = J -^ailj 



+ * 



(38) 



The fact that the coefficients a.^ do not depend on j relies 
on the translational invariance. All dimers are equal. 
But this is not necessary for our perturbation scheme to 
hold. On the contrary, we consider it one of the major 
advantages of the scheme presented here that it can be 
done in real space without knowing the form of the eigen 
states in advance. 

In 10 order the variable i in Eq. (pq) runs from -10 
to 10. If j is chosen appropriately, a chain segment with 
II dimers suffices to compute all the coefficients a^. For 
instance, for oio one has to take j = whereas for oq 
the right choice is j = 5. All sites that can be reached 
within 10 hops starting at j and ending at j -I- i must be 
contained in the cluster to avoid finite size effects. 

Of course, given translational invariance we know that 
spatial Fourier transform provides the eigen states |fc) = 
^2/_L^j exp{ikj)\j) characterized by their lattice mo- 
mentum k. The corresponding eigen energies read 



iu{k)^{k\H,ff\k)~Eo 



(39a) 



^ J -ao- Eo + J ■^2aj- cos{jk) . (39b) 

For dimerized chains the dispersion minimum, the 
triplet gap, is found at /c = 



A = u{k^O) 



(40) 



It is beyond the scope of the present paper to compute 
the effective Hamiltonian in the two-triplet sector. But 
an outlook on this issue is in order. Conventional per- 
turbation schemes have difficulties to compute proper- 
ties for two elementary particles because the structure of 
the eigen states is not known beforehand. In particular, 
bound states of two elementary particles are extremely 
difficult to obtain by conventional perturbative methods. 
But if the action of the effective Hamiltonian ( p5| ) on 
two triplets is calculated everything else can be deduced. 
First, one has to determine the coefficients A^-^^i^.j-^j^ de- 
fined by 



Hes\jij2) = X!^*i' 



>.\h,i2) , 



(41) 



where we assume that |ji, J2) has a triplet t^ on dimer j'l 
and a triplet t^^ on dimer J2- This is sufficient to com- 
pute triplets coupled to S'tot = I as was done successfully 
in Ref. ||7|]; the wave function ip{ji,J2) is antisymmetric 
under exchange ji <-> J2- Second, one has to diagonalize 
the matrix defined by the coefficients Ai-^i^.j-^j^. This 
can be carried out by standard Lanzcos algorithms. 

Spectral functions are also accessible by flow equation 
perturbation if the observable Q under study is unitarily 
transformed by the same transformation as the Hamilto- 



dQ{e) 
di 



[vie),Q{e)] 



(42) 



This is needed to know the matrix elements after the 
transformation. Preliminary studies showed that the 
treatment of Eq. (E3) is feasible within the perturbative 
approach. The necessary programmes are very similar to 
those for the Haniiltonian. Yet the treatment of observ- 
ables is more laborious than the one of the Hamiltonian 
itself. 



B. Computer Aided Evaluation 

Again C-l — h is the programming language of our choice. 
To encode the states of the dimerized chain it suffices to 
reserve two bits per dimer. For instance, four bytes can 
encode the state of a chain with 32 sites. The lowest 
bit represents site I, the second lowest represents site 2 
and so on. By applying the Ti in Eqs. ( ^ these states 
acquire polynomials in a as pref actors. Thus we choose 
the basic data elements to be objects of the class BDE 
as sketched in Fig. |^. 



class BDE: 



class POLY: 



unsigned int X; 
POLY *first: 
BDE *next; 



short int alphaexp; 
int prefact; 
POLY *next; 



FIG. 1. Sketch of the basic data element class 



The innermost part of the algorithm is the implemen- 
tation of the Ti matrices. Each matrix element is repre- 
sented by a block of C-l— I- code. These blocks allocate an 
appropriate number of basic data elements initialized by 
the current state of the chain. Then they modify those 
bits which represent the pair of adjacent dimers under 
study according to the rules of Table |. The Ti are im- 
plemented in five steps: 

(i) The number of sites L is chosen even. We introduce 
the loop variable 



pe{i,2, 



-} 



(43) 



such that p ~ 1 addresses dimer together with dimer 1, 
p = 2 addresses dimer 2 and dimer 1 and so on. Periodic 
boundary conditions are used, i.e. p = L/2 addresses 
dimer L/2 — 1 with dimer 0. 

(ii) From p a four byte long bit mask is constructed 
such that the bits refering to the dimers addressed 
by p are set to unity and zero otherwise. Thus 
p = 1 yields (1, 1, 1, 1,0, . . . ,0) and p = 2 yields 
(0,0,1,1,1,1,0, ...,0) and so on . 

(iii) The appropriate part x representing the two adjacent 
dimers is cut out by applying the logical AND operation: 
(bit mask) AND {X representing the chain state), 
(iv) The decimal value of x e {0, 1, ... , 15} is used to 
jump in a SWITCH(x)-CASE sequence to the appropri- 
ate CASE-block, encoding the matrix element of T,;. 
(v) The appropriate prefactor is multiplied to the allo- 
cated and modified data elements. 

Finally p is incremented by one and one goes back to (i) 
as long as p meets condition (B^) . 

Following these steps one obtains a sum of basic data 
elements describing an elementary chain state after the 
application of a T^. In general, the complete state after 
the application of a T^ is a linear combination of such 
elementary states encoded by the unsigned integers X. 
The linear combination is stored as a list. Since this 
list after its generation may contain identical states it is 
sorted according to increasing values of X by a quick- 
sort algorithm. Prefactors of identical states are added 
so that memory usage is reduced. 

To calculate products of Ti, the algorithm described 
above is applied repeatedly. If the chain state is a linear 
combination the algorithm is applied to each addend of 
this linear combination. 

To calculate the ground state energy (|3^) the expres- 
sion a''' J2\m\=k C{rn)T{m) from Eq. ( p5| ) must be applied 
to the triplet vacuum |0) for each order k. First the cur- 
rent index m and its coefficient C{rn) are read from an 
input file prepared previously, see Sect. II. B. Then T{m) 
is applied to |0) as described above. Since M{m) — 
T{m) reproduces |0) up to a prefactor which is a poly- 
nomial in a. This polynomial and C{m) are piped to 
an algebraic computer programme (Maple) to multiply 
them. This scheme is iterated with intermediate summa- 
tions by Maple till the final result is found. 



The result for the ground state energy per spin 
Hamiltonian (g) up to 10 order in A reads 



of the 






-2a) 



,A^- 



^89 

,48 



1^ 

311 



^3 

.4 



a 



-a - 



227 

257909 
^ 3456 



1307 2 .„ 3 
a^ - 42a^ 



93 o 45 ,\^5 /463 
4 2 y V 96 



13 27 

\ ^ 

16 4 

^5/463 



3 \_3 /13 27 



12 



■a 



879 5W7 
-——a-" A 



215995 2 
/ 414359 



159 .A -6 / 81557 

a A — 

2 y V 6912 

173579 . 14865 , 



432 
139801 



8477587 2 



V 12960 '^ 648 " "*" 12960 



152558 3 2774357 4 ^^^^ 5 4527 gA^s 



81 



-a 



1620 



-I- 



-a^ 



■a" 



/ 2354594813 7341879263 14053262981 
~ V 24883200 ^ 12441600 " ^ 6220800 ° 
1591335559 3 9560574943 4 8121212969 
345600 ^ 1555200 259200 

453741 (, 248391 7A-9 / 106469295871 

a^ H a^ A - 

64 32 y V 373248000 

82849717337 107584683283 2 89796462557 

o. -\ (y. -\ f 

46656000 15552000 5832000 

160938279937 4 57686123141 , 



46656000 
160938279937 

5832000 



iOOOZUUU 

57686123141 
972000 



143920286959 ^ 339171 , 336527 sA^io\ . ^ 
972000 "+^"-^^"^ ' (44) 



iA is used. 



where the shorthand A 

To calc ulat e the dispersion uj{k) the hopping elements 
Gi in Eq. ( |39b| ) have to be determined. The effect of T(to) 
with M{rri) = is to shift a triplet by at most \rn\ dimers. 
Hence it suffices to perform the calculation for a given a^ 
on an appropriate chain segment. The calculations are 
analogous to those for the ground state energy. Note that 
by Eq. (|3^ ) the general effect of a product T{m) will be 
a sum of states each containing one triplet on different 
dimers. For a given a^ the corresponding state has to 
be found in that list. Its polynomial prefactor yields a^. 
Results for the hopping elements up to 10 order in A 
are presented in appendix |Q. 

The calculations took 10h-30h strongly depending on 
the number of sites L; 100MB-500MB memory were used. 
The computations were done on a Sun Ultra workstation. 



IV. RESULTS AND COMPARISON TO OTHER 
METHODS 

So far all results refer to Hamiltonian (Q). In this sec- 
tion we prefer to present the results corresponding to 
Hamiltonian (0) by substituting according to Eq. (@). 
Since A is the expansion parameter it is substituted be- 
fore we manipulate the equations further, see below. The 



frustration a on the other hand is treated as a fixed pa- 
rameter throughout further manipulations. It is substi- 
tuted only at the very end. 



A. Ground State Energy 

Substituting J, A and a in Eq. ( p4|) leads to the dashed 
curve in Fig. called "plain series" . The solid line rep- 
resents a more sophisticated approach. Substituting J 
and A leads to an expression in 8 and a. In order to 
profit from the knowledge that the ground state energy 
is lowered by external dimerization as Ai?o oc —8^1'^ |15| ] 
(5 = a;2 is substituted. Thereby a is treated as fixed pa- 
rameter which does not change the exponents, here 4/3. 
Their values are protected by symm etry [|l6[ as long as 

This means 



logarithmic corrections are neglected 1 17 



)nenrs 



that this statement is true as long as a is below its crit- 
ical value ac = 0.241167(±5) [|l9|. Above this value the 
translational invariance is broken spontaneously and the 
ground state is dimerized so that the leading power be- 
comes linear Aii^o oc —S (see below). 
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- plain series 

- x=5" & Taylor in x, 
no linear x 




FIG. 2. Ground state energy per site e vs. dimerization 5 
without frustration (qq ~ 0) for Hamiltonian m). For details 
see main text. The inset shows an enlargement for small 5. 

The result for the ground state energy is expanded in 
a Taylor series in x up to order 10 about x = 1. This is 
the limit of strong dimerization for which our expansion 
holds. Adding the term Y ■ {1 — x)^^ introduces an un- 
known Y which can be fitted such that the linear term 
in X vanishes. In this way the leading order contribu- 
tion to the energy lowering agrees with the continuum 
prediction. The result of this procedure is depicted as 
solid curve in Fig. and is labeled accordingly. Both 
curves are compared to numeric results from Density Ma- 
trix Rcnormalization Group (DMRG) shown as symbols. 
The agreement is excellent. The accordance is better for 
the biased extrapolation based on the continuum theoret- 
ical power law as was to be expected. The extrapolated 
results can be trusted even quantitatively down to 1 or 
2%. 

In absence of frustration (a — ao ^ 0) the right hand 



side of Eq. ( p4| ) becomes a polynomial in A of degree 
10 with rational coefficients. Using a different perturba- 
tion method Barnes et al. [M calculated the same coef- 
ficients up to order 9, providing a good check. Further- 
more, Gelfand et al. |E1| also calculated a polynomial in 
A. Their method gives the coefficients as real numbers 
up to order 15. We could varify these numbers to the 
given precision in Ref. ||2l| up to order 10. 

One of the major motivations to study dimerized spin 
chains comes from spin-Peierls systems where spin and 
lattice degrees of freedom are coupled. In the adiabatic 
description of this phenomenon one adds the elastic en- 
ergy -y(5^ to the Hamiltonian (0). This term takes into 
account that it costs energy to modulate the magnetic 
couplings. To determine the equilibrium value of S the 
ground state energy is minimized by variation of 6 



= 



86 



KS . 



(45) 



So one has to know the derivative of eo which is given by 
the expectation value per site of the dimerization opera- 
tor i^DiM = iH{6) - H{0))/S with H from Eq. (|) 



{H^ 



dim; 



dS 



(46) 



In Fig. H we plot deo/dS as derived from our results for 
the ground state energy. Note that 3/8 is an upper bound 
for the expectation value of the dimerization since the 
dimerization is maximum if every second bond is occu- 
pied by a singlet (82^82^+1) = —3/4. This upper bound 
is excellently complied with by our results since we ex- 
pand around the limit of complete dimerization 6=1. 



0.40 




FIG. 3. Expectation value of the dimerization operator 
HuiM vs. dimerization S without frustration (ao = 0). 

Fig. |j show the ground state energy per site for 
ao = 0.241 and ao = 0.35. For supercritical frustra- 
tion we know that the leading term is linear in 6. In 
order to be able to describe well the crossover from (5 3 to 
6 behavior we use the substitution 6 = x^ and suppress 
linear and quadratic terms in the Taylor expansion by 
appropriate higher order terms in (1 — x). This proce- 
dure leads to a series in 6 comprising terms oc 6 and oc Ja 



as well as higher terms in 5^^^ . The agreement obtained 
in comparison to DMRG data is again excellent down to 
very low values of dimerization. In Fig. || the expectation 
values of the corresponding dimerization operators i?DiM 
are plotted. 
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FIG. 4. Ground state energy e vs. dimerization 5; frustra- 
tion ctQ = 0.241 (upper panel), cto — 0.35 (lower panel) 
. The insets show enlargements for small 6. 

We limited the range of 6 for finite frustration to the 
interval [0,0.2]. This is required by the nature of our 
perturbative expansion around isolated dimers. Fixing 
ao (not a !) at a finite value implies that the limit (5-^1 
is not the limit of isolated dimers. The point 6 = 1 cor- 
responds to a spin ladder with coupling 2 Jq on the rungs 
and agJo on the legs. Of course, our bare perturbative 
results pertain also to the ladder where they correspond 
to an expansion around the rung limit. This limit has 
already been investigated intensively |2^ so that we re- 
frain here from a comprehensive analysis. Note that the 
perturbation in the ladder is considerably simpler since 
iV = 1 in Eq. (||) whereas iV = 2 is treated here. In other 
words there is no creation or annihilation of two triplets 
but only shifts of triplets or creation or annihilation of 
one triplet. 

To understand the reason for the restricted applicabil- 
ity of the analysis presented here a simple comparison of 
couplings suffices. Using x = S"^ or x = S'^ is optimized 



to treat the case where the coupling Jo(l — 15) is the dom- 
inant perturbation. As a rule of thumb this is the case if 
2ao < 1 — (5. Otherwise one has to treat the frustrating 
coupling as the dominant perturbation. We restrict our 
present analysis to the regime where the 1 — S coupling is 
the dominant perturbation so that S may not be chosen 
too large. For ao ~ 0.35 the dimerization S should not 
exceed about 0.3. 
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FIG. 5. Expectation value of the dimerization operator 
-^DiM vs. dimerization S; frustration ao = 0.241 (upper 
panel), ao = 0.35 (lower panel). 



B. Energy Gap 

The first approach is again to substitute J, A and a 
in Eq. ( ^ according to Eq. (||). In the case ao = our 

result agrees to the 9 degree polynomial Barnes et al. 
[ p0| computed for the energy gap. For various values of 
ao the results are plotted as dashed curves labeled "plain 
series" in Figs. |[ ^ and ||. 

The solid curves were obtained by an extrapolation bi- 
ased by the continuum theory results A oc S'^^'^ for sub- 
critical frustration ||l^ and A — A|g^o oc S^'^ for super- 
critical frustration, see e.g. Ref. |23]. To this end, we 
substituted J and A in Eq. ( j^ ) according to Eq. (|^). 
Next (5 = X2 is replaced treating a as a fixed parameter. 



This result is re-expressed by an appropriate Pade ap- 
proximant about a; = 1, which corresponds to the dimer 
limit. As will be seen below the use of a Pade approx- 
imant instead of a series in x matters only for sizable 
frustration. The Pade approximant is chosen such that 
it takes the information about the 11 coefRcients of the 
expansion into account. In the final expression x — S'^ 
and a = -^ are inserted to obtain the data shown. The 
comparison to DMRG data [g^ shows that the biased 
extrapolation is extremely precise for most choices of pa- 
rameters. 
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FIG. 6. Energy gap A vs. dimerization (5 without frus- 
tration. The numbers in the square bracket stand for the 
polynomial degree in the numerator and the denominator of 
the Fade approximant, respectively. The inset shows an en- 
largement for small &. 
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FIG. 7. Energy gap A vs. dimerization 5 for critical frus- 
tration Qfo = 0.241. 

In order to check the reliability of our results in the 
most difficult case we extrapolate our results to zero 
dimerization. This means we employ the above proce- 
dure of a biased Pade approximant and set (5 = 0. The 
gap dependence on the frustration ao is depicted in Fig. y 
and compared to DMRG results 1 24 2^ . For comparison 
we include also an extrapolation of the series in a; = b'^l'^ . 
Rigorously, there is no gap below the critical value ac 
|P6| . The numerical value (oLf^ — 0.241 167(±5)) is taken 
from Ref. |19]. The wiggle at about (5 = 0.1 is a spuri- 



ous pole resulting from the Pade approximation which is 
not present in the series representation. The Pade ap- 
proximant, however, is better for larger frustration. The 
overall agreement is good but not excellent. Obviously, 
neither the Pade approximant nor the series in x are fit 
to describe the essential singularity at critical frustration. 
Yet, more sophisticated approximants which allow to ex- 
trapolated in several variables (here: A and ag) might 
render a more efficient analysis of the available expan- 
sion coefficients. 
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FIG. 8. Energy gap A vs. dimerization 5 for supercritical 
frustration ag = 0.35, oq = 0.5. 
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FIG. 9. Energy gap A vs. frustration a at 5 — 0. 
dashed vertical line indicates the value of Uc. 



The 



To complete the analysis of our expansion coefficients 
other approximants were also used. The common ap- 
proach to detect critical behaviour is to approximate the 
derived logarithm of the function under study by a Pade 
approximant (Dlog Pade). In this representation the po- 
sition Ac of a singularity (A — Xc)'^ is given by a pole and 
its residue defines the critical exponent 7. The results do 
not depend much on whether the Hamiltonian (0) or (0) 
is used. In fact, the results for Hamiltonian (Q) agree a 
bit better with the predictions Ac = 1 and 7 — 2/3 than 
those for Hamiltonian (Q). In absence of frustration we 
find Ac = 1.002 and 7 = 0.74. The value for the position 
is very encouraging; the result for the critical exponent, 
however, is a bit disappointing. The same observation 
was made in Ref. |Q. No significant change in the ex- 
ponent occurs if the position of the singularity Ac = 1 is 
pre-set. 

The reason for the difficulty to find the correct expo- 
nent is found in the logarithmic corrections. This can be 
seen in two ways. The first one is to go to the critical 
frustration where no logarithmic corrections are present. 
This proved useful in numerical analyses [p4| , too. In- 
deed, at ao = ttc we find 7 = 0.65 which agrees with 
2/3 within 3%. The second is to take the logarithmic 
corrections into account. Usually, they are computed for 
subcritical frustration and given in the form [E7 17|Jl^ 



(47) 



A{6)(xS^/y\\n{S/So)\^/^ 

where Sq is some unknown constant. Hence, we apply the 
Dlog Pade approximation to 

f{5) :- V-M'5) + M«o)A , (48) 

where h{ao) is some constant depending on the frustra- 
tion. We look at the biased approximant /'// = -P(^^)/^ 
about S^ — 1 where -P(5^) is a polynomial of order 9 in 
(5^ which takes the obvious symmetry 6 ^^ —6 oi Hamil- 
tonian (|^) into account. To be able to expand about 
isolated dimers the replacement ag —> ao{l — S'^) is car- 
ried out. The value P{0) then is a direct estimate for the 
critical exponent 7. Assuming that the coefficients of P 
decrease quickly on increasing order once the logarithmic 
correction is properly taken into account we fix the con- 
stant h{ao) such that the coefficient of order 9 vanishes. 
In this way we find 7 = 0.68 and h = 5.65 at zero frustra- 
tion. At critical frustration the values are 0.56 and 2.50, 
respectively. 

The good agreement without frustration indicates that 
logarithmic corrections are indeed the reason for the dif- 
ficulty to determine the correct exponents from the per- 
turbative data. The unsatisfactory agreement at critical 
frustration tells us that the way h is determined is not 
optimum. It would be very helpful if a prediction for the 
gap as function of S and a — ac existed which comprised 
also the regime of supercritical frustration. 



C. Dispersions 

The one magnon dispersion uj(k) for t he an tiferromag- 
netic Heisenberg chain is given by Eq. (p9b|). The hop- 
ping elements a^ are given in appendix ICl. For ag = 
they can be compared till order 5 with those computed 
in Ref. IpO]. Full agreement is found. 

To deduce the dispersion relations from the bare coef- 
ficients several approaches will be presented. T he di rect 
approach is again to substitute J, A and a in Eq. ( |39h| ) ac- 
cording to Eq. (P). By construction, the resulting curves 
reproduce the same energy gaps at k — as the gap re- 
sults labeled "plain series" in the preceding section. We 
learned, however, in the preceding section that a biased 
extrapolation is very useful to approximate the gaps. So 
the question arises how the good biased extrapolations 
can be used for the description of dispersion relations. 

Motivated by the behaviour of Lorentz-invariant sys- 
tems where the dispersion passes from uj = vgk (wg spin 
wave velocity) to lj ex \/A'^ + (wsfc)^ when a gap opens 
we use th e fo llowing procedure. We substitute J and A 
into Eq. (B9q) leading to the plain series result we shall 
refer to as Wpiain(fc; a)- Then the difference 



■^diff 



(fc) 



,(fc;a) -Wpiai„(fc = 0;a) 



(49) 



is expanded as Taylor series about (5 = 1 up to 10 order 
while a is treated as fixed parameter. Only then a is 
substituted according to Eq. (^. Finally, the dispersion 
is computed by the quadratic mean 



Lo{k) - ^A^{ao,6)+Lul^{k) 



(50) 



For 6 — ao — the result is shown in Fig. (|T^) where it is 
also compared to the plain series result to illustrate the 
effect of the quadratic mean. The curve is compared to 
the rigorous result known from Bethe ansatz ||28| . The 
quadratic mean matters only for low energies. Figs. O, 
[1^ and |lj show the corresponding results for various S 
and a values for relatively low frustration in the subcrit- 
ical regime or close to it. As soon as there is some dimer- 
ization, for instance 5 > 0.05 in Fig. O, the difference 
between the plain series (not shown) and the quadratic 
mean (pO) is no longer discernible. 
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FIG. 10. Dispersion w vs. wave vector k without frustra- 
tion and dimerization.The curve Wpiain depicts the plain series 
result. 
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FIG. 11. Dispersion u vs. wave vector fc at ao = and 
various dimerizations. 
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FIG. 12. Dispersion u vs. wave vector fc at ao = 0.241 and 
various dimerizations. 



1.50 

1.25 
^ 1.00 
O 0.75 

O 

3 0.50 

0.25 

0.00 

0.0 0.2 



5=0 


8=0.05 


8=0.1 ^--^J^S^JZZ 


8=0.2 ^--'''^-<5?^^^^''" 








_--'"' ^^^.^'^ 


.^y""''^ ■'^ -y' 












""""^ ^^ y 


."' /' 


,^'' ^ 


.--'' y 


y 


y 


y 


y 


1 . 1 . 1 . 1 ■ 



0.4 0.6 

k/TT 



1.0 



FIG. 13. Dispersion ci; vs. wave vector fc at ao = 0.35 and 
various dimerizations. 

Entering the regime of large frustration (ag > 0.4) the 
calculations have to be modified slightly. In Fig. 11^ we 
compare two different approaches to calculate the disper- 



sion Lo(k'). The dashed line labeled "wpiain(fc) & Taylor" 
is obtained by expanding a;piain(fc;a) in a Taylor series 
about b — \ with a as fixed parameter. Finally, a is 
substituted according to Eq. (^) for 5 = 0, an = 0.5. The 
resulting curve is far off for low wave vectors. To get good 
values for the energy gap it is necessary to use the square 
root representation (pO[). On the other hand, the value 
at fc = TT is unity which is the exact result [£9|,Q . If the 
square root representation (50) is used the dashed-dotted 
line labeled "sqrt & Taylor" in Fig. |lj is obtained. The 
energy gap is nicely reproduced BQl but the agreement 
at the dispersion maximum deteriorates considerably. 

So one would like to have a representation at hand 
which combines the advantages of the two approaches. 
A convincing interpolation is achieved by replacing the 
Taylor series by the corresponding [5,5] Fade approxi- 
mant labeled "sqrt & Fade" in Fig. n4l Fig. 15 presents 



in Fig. y_4[ Fig. 
the results of the latter approach for finite b. For com- 
pleteness, we mention that in the ap ^ 0.4 regime the 
Fade representation works as well as the Taylor expan- 
sion. For instance, the Fade representation reproduces 
the results presented in Fig. O within 1%. This means 
that for ao < 0.4 there is no need to use the more tedious 
Fade representation. But its use for larger frustrations 
does not imply an inconsistency of our general approach. 
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FIG. 14. Dispersion uj vs. wave vector A; at ao = 0.5 and 
(5 = 0. The results of three different ways to deduce the 
dispersion relation from the bare coefficients are depicted. 
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FIG. 15. Dispersion a; vs. wave vector A; at ao = 0.5 and 
various dimerizations. The approach used is the one labeled 

The observation that the dispersion relation for strong 
frustration behaves at fc « and at A: w tt quaUtatively 
differently can be understood on physical grounds. For 
this purpose let us consider the Majumdar- Ghosh model 
at ao = 1/2 without dimerization. Its low-lying excita- 
tions are asymptotically free S = 1/2 spinous |pl|- p3| , p0[ . 
This implies that in the vicinity of the dispersion mini- 
mum the dispersion relation for S ^ does not represent 
the branch of a well-defined magnon excitation. But it is 
the lower band edge of a two-spinon continuum. In this 
respect it is similar to the situation of subcritical frustra- 
tion. So we have to use the corresponding appropriate 
extrapolation to obtain reliable results. 

In the vicinity of the dispersion maximum, however, 
the spinon interaction dominates of their kinetic energy 
and binding occurs even without dimerization. For in- 
stance, the exact triplet state p9[ is a tightly bound 
two-spinon state. A variational analysis shows that this 
bound triplet exists below the two-spinon continuum in 
a finite interval around k = n |pl| , p3 30 1. Hence, it 
is fully comprehensible that the perturbative approach 
which starts from local triplets works more easily around 
k = TT whereas improved extrapolation is necessary 
around k = 0. 



V. CONCLUSIONS 

In summary, we presented a perturbative scheme which 
relies on a suitably chosen unitary transformation. The 
scheme works for an unperturbed equidistant spectrum 
which can be labeled by the number of energy quanta. 
The perturbation term changes the number of energy 
quanta at maximum by a finite number N. Thereby, 
we generalized the approach by Stein {N = 1) ||^ to gen- 
eral N. The unitary transformation is carried out by flow 
equations Q. 

By the transformation a systematic mapping of the 
original problem to an effective Hamiltonian is achieved 
which conserves the number of energy quanta. Thus 
Hilbert space sectors with different number of energy 
quanta are separated. Our scheme will be particularly 
useful where ordinary perturbation theory is hampered 
by the fact that the structure of the perturbed states 
is not known, for instance, systems without translation 
invariance or two-particle problems. The possibility to 
obtain information in symbolic form, i.e. as polynomi- 
als, distinguishes our approach from the multiprecision 
method by Barnes et al. [E0[ and other conventional im- 
plementations [ pl[ . 

The realization of the perturbative scheme comprises 
two distinct steps. The first is still general, the second 
specific to the model. For the first step we provided the 
necessary coefficient up to order 10 for N = 2. In written 



form they are included till order 6 in this publication. 
The other coefficients shall be provided electronically. 

The second step is illustrated by dimerized frustrated 
S = 1/2 chains. The limit of isolated dimers has an 
equidistant spectrum. The ground state is the product 
of singlets on the dimers. The energy quanta are triplets 
on the dimers. The perturbing weak couplings between 
the dimers may create/annihilate at most two triplets so 
that N — 2 holds. After the transformation the num- 
ber of triplets is conserved. We presented results for the 
ground state (zero triplets) and the magnon dispersion 
(one triplet). Thereby we demonstrated the validity and 
applicability of the scheme proposed. 

The results for the spin chains are given as polynomials 
in the frustration parameter a. So they are easy to use for 
anybody who wants to analyze data by appropriate fits. 
No new calculations are necessary. The fits can be carried 
out instantly. Based on the results for the dimerization 
operator Hy) the dependence of the dimerization 6 on the 
elastic spring constant K is quickly accessible. 

As we demonstrated our results are reliable down to 
about 6% without additional information, i.e. using the 
plain series. With additional information as the criti- 
cal exponent, for instance, the results can reliably be 
used down to about 2%, in some cases even less. More- 
over, starting from the exact coefficients more elaborate 
schemes like differential approximants in two variables 
become possible P4[] . 

As an outlook we like to point out that our approach 
can also be used to compute dispersion relations in two 
or higher dimensional dimerized spin systems as demon- 
strated recently for (VO)2P207 Q. Investigations for 
CuGeOa are in preparation. They will improve consid- 
erably the third order analysis ||35| . Another fascinating 
field concerns computations in the two-magnon sector. 
The attractive interaction of two magnons pq |, for in- 
stance, leads to bound states. For (VO)2P207, which 
is characterized by a relatively large dimerization, such 
results were found in a fourth order calculation. If the 
unitary transformation is applied to observables like, for 
instance, the Raman operator spectral functions are also 
within reach. 
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It is shown that the choice ( |12D for the infinitesi- 
mal generator achieves block diagonaUty for the eS^ective 
Hamiltonian H{£ = cxo). The proof follows the lines of 
Mielkes proof for band matrices WM. Let {|i^i)} be the 
eigen state basis of Hq and define 






Inspecting Eqs. (^;|0 12) closely one realizes that the 
choice for the infinitesimal generator is equivalent to 

r,.,{i)=sgn{hl-h%)h,,{e). (Al) 

Inserting this expression in flow equation (|g) yields 

+ J2 (sgn(/»° - ht,) + sgn{h% - hl,))\h,k\' . (A2) 

Assume without loss of generality that the eigen states 
\i^i) are labeled such that /i"^, > h% ii k > i. Let us 
consider the sum of the first r diagonal elements of H{£) 



d 



-J2hu^2j2J2'Sn{hl~hl,)\h,,\ 



I— 1 



(A3) 



i—l k>r 



The right side of Eq. (A3) is non-positive. Thus the sum 
on the left hand side is a continuous monotonically de- 
creasing function with £. If we know that it is bounded 
from below we conclude that the sum converges whence 
the vanishing of its derivative for ^ — > 00 ensues imme- 
diately. Hence we need beyond the conditions (i) and 
(ii) in the Introduction the boundedness for the whole 
Hamiltonian. If this is not given our choice for 77 might 
be problematic whereas Wegner's choice still works (for 
an example, see Ref. ]36|). 

As we are, however, interested in deriving a perturba- 
tion expansion order by order we can assume the whole 
Hamiltonian to be bounded from below without loss of 
generality. For any finite order of the expansion it is suf- 
ficient to consider a finite cluster supposing some short 
range interaction. Then the Hamiltonian is a finite di- 
mensional matrix and is certainly bounded from below. 
This is true in particular if we stay on the abstract level 
as in Eq. (0). The generalized variational principle im- 
plies that 



H'^-^II' 



(A4) 



holds for all £ where the cui are the eigen values of H in 
ascending order. Note that the eigen values are invariant 
under the unitary transformation. 
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The vanishing of the derivative of J2l=i ^a ^'^^ arbi- 
trary r impUes eventually 

\imsgn{h%^hl,)\h,k\'=0. (A5) 



From this equation follows that either the cigen states 
\i>i) and Ivk) are degenerate, i.e. they belong to the same 
block, or hik{i = oo) = 0, i.e. matrix elements link- 
ing different blocks vanish. Hence block diagonality is 
achieved and the number of energy quanta given by Hq 
is conserved 

\im[Ho,H{£)]=0 . (A6) 

This concludes the formal proof. 

To restrict the argument at one stage to finite clusters 
does not constitute a real restriction for the series expan- 
sion in any finite order. The linked cluster theorem tells 
us that any finite order can be found from an appropriate 
finite cluster. 

In practice, the important issues is whether the physics 
remains the same on variation of the expansion parame- 
ter A. In our example of dimerized chains the expansion 
makes sense as long as the gap does not close. In other 
words, those situations are accessible which can be linked 
continuously to the dimer limit by gapped systems. Gap- 
less situations can only be described if they are the limit 
of gapped systems linked to the dimer limit. 
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APPENDIX B: COEFFICIENTS OF THE EFFECTIVE HAMILTONIAN 



In the following tables the coefficients C{rn) of the effective Hamiltonian 
Order 7 to 10 will be available electronically. 



5[) are given up to sixth order inclusively. 



m 


C{m) 


m = 1 





1 


\rn\ = 2 


1-1 


1 


2-2 


1/2 


TO = 3 


01-1 


-1/2 


02-2 


-1/8 


10-1 


1 


11-2 


1/2 


1-21 


-1 


20-2 


1/4 


TO — 4 


001-1 


1/4 


002-2 


1/32 


010-1 


-1 


011-2 


-3/8 


01-21 


1/4 


01-10 


1/2 


020-2 


-1/8 


02-20 


1/16 


02-1-1 


-3/8 


100-1 


1 


101-2 


1/2 


10-21 


-1 


110-2 


1/4 


11-1-1 


1/2 


12-2-1 


1/3 


12-1-2 


1/6 


1-22-1 


-1 


1-2-12 


1/2 


1-11-1 


-1 


1-12-2 


-3/8 


1-1-22 


1/8 


200-2 


1/8 


21-1-2 


1/12 


22-2-2 


1/16 


2-22-2 


-1/8 


2-11-2 


1/4 


|to| = 5 


0001-1 


-1/8 


0002-2 


-1/128 


0010-1 


3/4 


0011-2 


7/32 


001-21 


-5/16 


001-10 


-3/8 


0020-2 


3/64 


002-20 


-3/128 


002-1-1 


7/32 


0100-1 


-3/2 


0101-2 


-5/8 



TO 


C(to) 


TO = 5 


010-21 


7/8 


010-10 


3/2 


0110-2 


-1/4 


011-20 


7/16 


011-1-1 


-5/8 


012-2-1 


-7/18 


012-1-2 


-11/72 


01-201 


1/8 


01-210 


-3/8 


01-22-1 


1/2 


01-2-12 


-1/24 


01-11-1 


9/8 


01-12-2 


35/96 


01-1-22 


-19/96 


01-1-11 


-3/8 


0200-2 


-3/32 


020-20 


3/32 


020-1-1 


-1/4 


021-2-1 


-11/72 


021-1-2 


-1/18 


022-2-2 


-5/128 


02-21-1 


7/24 


02-22-2 


9/128 


02-2-22 


-3/128 


02-2-11 


-1/12 


02-10-1 


-5/8 


02-11-2 


-1/4 


02-1-21 


7/24 


1000-1 


1 


1001-2 


1/2 


100-21 


-1 


1010-2 


1/4 


101-1-1 


1/2 


102-2-1 


1/3 


102-1-2 


1/6 


10-201 


1 


10-22-1 


-1 


10-2-12 


1/2 


10-11-1 


-3/2 


10-12-2 


-5/8 


10-1-22 


3/8 


10-1-11 


1/2 


1100-2 


1/8 


110-1-1 


1/4 


111-2-1 


1/6 


111-1-2 


1/12 


112-2-2 


1/16 


11-21-1 


-5/8 


11-22-2 


-1/4 


11-2-22 


1/8 



TO 


C(to) 


TO = 5 


11-2-11 


1/8 


11-11-2 


1/4 


11-1-21 


-1/2 


120-2-1 


1/9 


120-1-2 


1/18 


121-2-2 


1/24 


12-21-2 


1/6 


12-2-21 


-1/3 


12-10-2 


1/12 


1-202-1 


1 


1-20-12 


-1/2 


1-211-1 


1/4 


1-212-2 





1-21-22 


1/4 


1-21-11 


3/4 


1-221-2 


-1/2 


1-22-21 


1 


1-2-212 


-1/6 


1-2-102 


-1/4 


1-101-1 


3/4 


1-102-2 


7/32 


1-10-22 


-3/32 


1-10-11 


-1/4 


1-111-2 


-5/8 


1-120-2 


-1/4 


1-1-202 





1-1-1-12 


1/8 


2000-2 


1/16 


201-1-2 


1/24 


202-2-2 


1/32 


20-22-2 


-3/32 


20-2-22 


1/32 


20-11-2 


1/8 


210-1-2 


1/36 


211-2-2 


1/48 


21-21-2 


1/12 


220-2-2 


1/64 


2-202-2 


3/64 


2-20-22 


-1/64 


2-211-2 


-1/4 


2-2-1-12 


1/8 


2-101-2 


1/4 


2-1-2-12 


1/4 


TO = 6 


00001-1 


1/16 


00002-2 


1/512 


00010-1 


-1/2 


00011-2 


-15/128 


0001-21 


5/64 


0001-10 


1/4 



TO 


C{m) 


TO = 6 


00020-2 


-1/64 


0002-20 


1/128 


0002-1-1 


-15/128 


00100-1 


3/2 


00101-2 


17/32 


0010-21 


-3/4 


0010-10 


-3/2 


00110-2 


11/64 


0011-20 


-45/128 


0011-1-1 


17/32 


0012-2-1 


17/54 


0012-1-2 


85/864 


001-201 


1/2 


001-210 


1/64 


001-22-1 


-13/24 


001-2-12 


53/288 


001-100 


3/8 


001-11-1 


-29/32 


001-12-2 


-301/1152 


001-1-22 


179/1152 


001-1-11 


7/32 


00200-2 


3/64 


0020-20 


-3/64 


0020-1-1 


11/64 


0021-2-1 


85/864 


0021-1-2 


43/1728 


0022-2-2 


17/1024 


002-200 


3/256 


002-21-1 


-197/1152 


002-22-2 


-29/1024 


002-2-22 


7/1024 


002-2-11 


-23/1152 


002-10-1 


17/32 


002-11-2 


11/64 


002-1-21 


-65/288 


002-1-10 


-45/128 


01000-1 


_2 


01001-2 


-7/8 


0100-21 


23/16 


0100-10 


3 


01010-2 


-3/8 


0101-20 


17/16 


0101-1-1 


-7/8 


0102-2-1 


-5/9 


0102-1-2 


-17/72 


010-201 


-5/8 


010-210 


-3/4 


010-22-1 


9/8 


010-2-12 


-7/18 


010-11-1 


39/16 



TO 


C(to) 


TO = 6 


010-12-2 


269/288 


010-1-22 


-181/288 


010-1-11 


-17/16 


01100-2 


-5/32 


0110-20 


11/32 


0110-1-1 


-3/8 


0111-2-1 


-17/72 


0111-1-2 


-7/72 


0112-2-2 


-9/128 


011-21-1 


5/6 


011-22-2 


113/384 


011-2-22 


-67/384 


011-2-11 


-7/24 


011-10-1 


-7/8 


011-11-2 


-3/8 


011-1-21 


7/12 


011-1-10 


17/16 


0120-2-1 


-4/27 


0120-1-2 


-13/216 


0121-2-2 


-25/576 


012-20-1 


-5/9 


012-21-2 


-17/72 


012-2-21 


103/288 


012-2-10 


17/27 


012-10-2 


-7/72 


012-1-20 


85/432 


012-1-1-1 


-17/72 


01-2001 


-9/16 


01-202-1 





01-20-12 


-1/6 


01-211-1 


-2/3 


01-212-2 


-15/64 


01-21-22 


37/192 


01-21-11 


13/24 


01-220-1 


7/8 


01-221-2 


1/3 


01-22-21 


-5/16 


01-22-10 


-11/12 


01-2-212 


-1/144 


01-2-221 


-25/288 


01-2-102 


-1/36 


01-2-111 


-1/6 


01-2-120 


31/144 


01-101-1 


-1 


01-102-2 


-329/1152 


01-10-22 


53/384 


01-10-11 


3/16 


01-110-1 


33/16 


01-111-2 


73/96 


01-11-21 


-49/48 



TO 


C(to) 


TO = 6 


01-11-10 


-27/16 


01-120-2 


151/576 


01-12-20 


-149/384 


01-12-1-1 


73/96 


01-1-202 


47/576 


01-1-211 


17/96 


01-1-220 


113/1152 


01-1-101 


7/16 


01-1-12-1 


-31/48 


01-1-1-12 


17/96 


02000-2 


-1/16 


0200-20 


3/32 


0200-1-1 


-5/32 


0201-2-1 


-7/72 


0201-1-2 


-11/288 


0202-2-2 


-7/256 


020-21-1 


73/288 


020-22-2 


39/512 


020-2-22 


-17/512 


020-2-11 


-17/288 


020-10-1 


-3/8 


020-11-2 


-5/32 


020-1-21 


2/9 


0210-2-1 


-13/216 


0210-1-2 


-5/216 


0211-2-2 


-19/1152 


021-20-1 


-17/72 


021-21-2 


-7/72 


021-2-21 


19/144 


021-10-2 


-11/288 


021-1-20 


43/864 


021-1-1-1 


-7/72 


0220-2-2 


-3/256 


022-20-2 


-7/256 


022-2-20 


17/512 


022-2-1-1 


-9/128 


022-1-2-1 


-25/576 


022-1-1-2 


-19/1152 


02-201-1 


-211/1152 


02-202-2 


-1/32 


02-20-22 


3/512 


02-20-11 


-7/384 


02-210-1 


95/144 


02-211-2 


85/384 


02-21-21 


-53/192 


02-220-2 


33/512 


02-22-20 


-27/512 


02-22-1-1 


85/384 


02-2-202 


7/512 


02-2-211 


23/384 
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TO 


C{m) 


TO 


C(to) 


TO 


C(to) 


TO 


C(to) 


TO 


C(to) 


TO 


C(to) 


TO = 6 


TO = 6 


TO = 6 


TO = 6 


TO — 6 


TO = 6 


02-2-101 


31/144 


10-212-2 


3/16 


11-2-12-1 


1/4 


1-20-212 


1/6 


1-12-2-11 


-3/16 


21-1-2-22 


5/288 


02-2-12-1 


-11/64 


10-21-22 


1/16 


11-101-2 


1/4 


1-20-102 


1/4 


1-12-11-2 


-3/8 


21-1-11-2 


1/24 


02-2-1-12 


23/384 


10-21-11 


3/8 


11-10-21 


-1/2 


1-2101-1 


-7/16 


1-1-2002 


-1/32 


2200-2-2 


1/256 


02-100-1 


-7/8 


10-220-1 


-1 


11-110-2 


1/8 


1-2102-2 


-5/32 


1-1-21-12 


-1/8 


221-1-2-2 


1/320 


02-101-2 


-3/8 


10-221-2 


-1/2 


11-11-1-1 


1/4 


1-210-22 


7/32 


1-1-221-1 


-11/48 


222-2-2-2 


1/384 


02-10-21 


7/12 


10-22-21 


1 


11-12-2-1 


1/6 


1-210-11 


11/16 


1-1-222-2 


-25/384 


22-22-2-2 


1/128 


02-110-2 


-5/32 


10-2-212 


-1/6 


11-12-1-2 


1/12 


1-2111-2 


1/4 


1-1-22-22 


19/384 


22-2-22-2 


-7/256 


02-11-20 


11/32 


10-2-221 


-1/3 


11-1-22-1 


-1/2 


1-2120-2 


1/16 


1-1-2-222 


-1/128 


22-2-2-22 


3/256 


02-11-1-1 


-3/8 


10-2-102 


-1/4 


11-1-2-12 


1/4 


1-21-202 


-3/16 


1-1-2-112 


-1/72 


22-2-11-2 


1/32 


02-12-2-1 


-17/72 


10-2-111 


-1/2 


11-1-11-1 


-7/8 


1-21-12-1 


1/2 


1-1-10-12 


-3/8 


22-11-2-2 


1/192 


02-12-1-2 


-7/72 


10-101-1 


3/2 


11-1-12-2 


-3/8 


1-21-1-12 


-1/2 


1-1-111-1 


-1/2 


22-1-21-2 


1/48 


02-1-201 


-5/36 


10-102-2 


17/32 


11-1-1-22 


1/4 


1-2201-2 


-1/2 


1-1-112-2 


-19/96 


2-2002-2 


-1/64 


02-1-22-1 


5/12 


10-10-22 


-9/32 


11-1-1-11 


3/8 


1-2210-2 


-1/4 


1-1-11-22 


23/96 


2-2011-2 


11/64 


02-1-2-12 


-11/96 


10-10-11 


-1/2 


1200-2-1 


1/27 


1-222-1-2 


-1/6 


1-1-121-2 


5/24 


2-20-1-12 


5/64 


02-1-11-1 


5/6 


10-110-1 


_2 


1200-1-2 


1/54 


1-22-22-1 


1 


1-1-1-212 


-5/72 


2-2101-2 


-3/8 


02-1-12-2 113/384| 


10-111-2 


-7/8 


1201-2-2 


1/72 


1-22-2-12 


-1/2 


1-1-1-102 


-1/8 


2-21-2-12 


-1/12 


02-1-1-22 


-67/384 


10-11-21 


5/4 


120-21-2 


1/18 


1-22-11-1 


1 


20000-2 


1/32 


2-21-12-2 


49/192 


02-1-1-11 


-7/24 


10-120-2 


-3/8 


120-2-21 


-1/9 


1-22-12-2 


3/8 


2001-1-2 


1/48 


2-21-1-22 


-9/64 


10000-1 


1 


10-12-1-1 


-7/8 


120-10-2 


1/36 


1-22-1-22 


-1/8 


2002-2-2 


1/64 


2-22-22-2 


1/16 


10001-2 


1/2 


10-1-202 


-1/8 


1210-2-2 


1/96 


1-22-1-11 





200-22-2 


-1/16 


2-22-2-22 


-3/128 


1000-21 


-1 


10-1-211 


-1/8 


121-20-2 


1/48 


1-2-2012 


1/18 


200-2-22 


1/32 


2-22-11-2 


-5/32 


10010-2 


1/4 


10-1-12-1 


3/4 


121-1-2-1 


1/36 


1-2-2102 


1/12 


200-11-2 


1/16 


2-2-21-12 


-3/32 


1001-1-1 


1/2 


10-1-1-12 


-1/8 


121-1-1-2 


1/72 


1-2-222-1 


-1/3 


2010-1-2 


1/72 


2-2-222-2 


-1/64 


1002-2-1 


1/3 


11000-2 


1/16 


122-2-2-1 


1/45 


1-2-22-12 


1/6 


2011-2-2 


1/96 


2-2-10-12 


-1/4 


1002-1-2 


1/6 


1100-1-1 


1/8 


122-2-1-2 


1/90 


1-2-2-122 


1/24 


201-21-2 


1/24 


2-2-112-2 


-17/192 


100-201 


1 


1101-2-1 


1/12 


122-1-2-2 


1/120 


1-2-1002 


1/8 


201-10-2 


1/48 


2-2-121-2 


1/24 


100-22-1 


-1 


1101-1-2 


1/24 


12-201-2 


1/6 


1-2-112-1 


-1/2 


2020-2-2 


1/128 


2-1001-2 


1/4 


100-2-12 


1/2 


1102-2-2 


1/32 


12-20-21 


-1/3 


1-2-11-12 


1/4 


202-20-2 


1/64 


2-10-2-12 


1/4 


100-11-1 


-2 


110-21-1 


-3/8 


12-210-2 


1/12 


1-2-121-1 


1/24 


202-1-1-2 


1/96 


2-11-11-2 


1/8 


100-12-2 


-7/8 


110-22-2 


-5/32 


12-22-2-1 


1/9 


1-2-122-2 


1/12 


20-202-2 


3/64 


2-12-21-2 


1/12 


100-1-22 


5/8 


110-2-22 


3/32 


12-22-1-2 


1/18 


1-2-12-22 


-5/24 


20-20-22 


-1/64 


2-1-221-2 


-1/4 


100-1-11 
10100-2 


1 

1/8 


110-2-11 
110-11-2 


1/8 
1/8 


12-2-22-1 
12-2-2-12 


-1/3 
1/6 


1-2-12-11 
1-2-1-222 


-13/24 
1/16 


20-211-2 
20-220-2 


-5/32 
-1/16 


2-1-111-2 


-3/8 






1010-1-1 


1/4 


110-1-21 


-1/4 


12-2-11-1 


-5/9 


1-2-1-112 


1/12 


20-2-1-12 


1/32 






1011-2-1 


1/6 


1110-2-1 


1/18 


12-2-12-2 


-17/72 


1-1001-1 


-1/2 


20-101-2 


1/8 






1011-1-2 


1/12 


1110-1-2 


1/36 


12-2-1-22 


11/72 


1-1002-2 


-15/128 


20-110-2 


1/16 






1012-2-2 


1/16 


1111-2-2 


1/48 


12-2-1-11 


2/9 


1-100-22 


5/128 


20-12-1-2 


1/24 






101-21-1 


-7/8 


111-21-2 


1/12 


12-100-2 


1/24 


1-1011-2 


17/32 


20-1-2-12 


1/8 






101-22-2 


-3/8 


111-2-21 


-1/6 


12-11-2-1 


1/18 


1-1020-2 


11/64 


20-1-12-2 


-5/32 






101-2-22 


1/4 


111-10-2 


1/24 


12-11-1-2 


1/36 


1-10-202 


3/64 


20-1-1-22 


3/32 






101-2-11 


3/8 


111-1-1-1 


1/12 


12-12-2-2 


1/48 


1-10-1-12 


5/32 


2100-1-2 


1/108 






101-10-1 


1/2 


1120-2-2 


1/64 


12-1-21-1 


-17/72 


1-1101-2 


-7/8 


2101-2-2 


1/144 






101-11-2 


1/4 


112-20-2 


1/32 


12-1-22-2 


-7/72 


1-1110-2 


-3/8 


210-21-2 


1/36 






101-1-21 


-1/2 


112-2-1-1 


1/16 


12-1-2-22 


1/18 


1-112-1-2 


-17/72 


2110-2-2 


1/192 






1020-2-1 


1/9 


112-1-2-1 


1/24 


12-1-2-11 


5/72 


1-11-2-12 


-7/24 


211-1-1-2 


1/144 






1020-1-2 


1/18 


112-1-1-2 


1/48 


12-1-11-2 


1/12 


1-11-11-1 


2 


212-2-1-2 


1/180 






1021-2-2 


1/24 


11-201-1 


17/32 


12-1-1-21 


-1/6 


1-11-12-2 


71/96 


212-1-2-2 


1/240 






102-20-1 


1/3 


11-202-2 


11/64 


1-2002-1 


-1 


1-11-1-22 


-43/96 


21-201-2 


1/12 






102-21-2 


1/6 


11-20-22 


-5/64 


1-200-12 


1/2 


1-11-1-11 


-3/4 


21-22-1-2 


1/36 






102-2-21 


-1/3 


11-20-11 


-5/32 


1-2011-1 


3/8 


1-1200-2 


-5/32 


21-2-2-12 


1/12 






102-10-2 


1/12 


11-211-2 


-3/8 


1-2012-2 


5/16 


1-121-1-2 


-7/72 


21-2-12-2 


-7/72 






102-1-1-1 


1/6 


11-21-21 


1/2 


1-201-22 


-9/16 


1-122-2-2 


-9/128 


21-2-1-22 


1/18 






10-202-1 


1 


11-220-2 


-5/32 


1-201-11 


-11/8 


1-12-21-1 


31/48 


21-11-1-2 


1/72 






10-20-12 


-1/2 


11-22-1-1 


-3/8 


1-2021-2 


1/2 


1-12-22-2 


83/384 


21-12-2-2 


1/96 






10-211-1 


5/8 


11-2-202 


-1/32 


1-202-21 


-1 


1-12-2-22 


-41/384 


21-1-22-2 


-11/288 
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APPENDIX C: EFFECTIVE HOPPING ELEMENTS FOR THE SPIN CHAIN 

In the following the effective hopping elements as they appear in Eq. ( p^ ) are given up to order 6. The effective 
hopping elements up to order 10 will be provided electronically. We substituted a = 1 — 2a and A = ^A. 

ao - So = 1 - 1 4 - Sa^ ] A^ - (8 - 8a - 6a^ + 3aAx^ - [2 - 24a + Sa^ + 8a^ + — a'' ] A* + (56 - 82a - 22a^ 



rr -s o^-4 r.^ .,x" ,'367 7328 22976 o 6442 . 28895_4 „^ , ,„_fi\^6 
SSa'' - 39a^ + 20a'' ) A +\ — —a + -^^a^ + ^^" ^J~'^ + ^^^'^ ~ ^^'^ ' ^ 

ai = -2aA - 4A^ - ( 8 - 8a - 25^ j A% ( 4 + 20a - 24^^ + lOo^ - 5a'' jA^ + [92 - —-a - i— a^ + 1520^ 

.„-A 13 5\^5 /532 11906 II960 , 1648 , 41357 , „, , r.^\r& 

470* H a^ I A + a H a^ + a^ a* + 85a^ + 6a^ | A 

\ 3 9 9 3 54 

o— z ; o i\— 3 / 9 c! 1 A— 4 / 272 220 , II50 o / 



aa = -a-'A - I 4a^ - 2a-' 1 A + I 6 - 4a - 23a^ + 140-^ - -cT 1 A + ( 36 - -^a - -^W + —^TT - 9cr 
f3 A-5 /I07 1630 1126 o 5102 . I3205_4 ^„ , ,,_6ix" 

aa = -a^A^ - {'^a' + 4a^ - 2a^>^ - f'-^a + 20a^ - ^^3 _ ,,^ _ ,-,\j^ _ /58 ^ 104_ _ 224_, __ ^^_, 
103_4 57 , 81_e\-6 

5^T^ _ ^40^3 ^ «^ _ ,^5\t^ ^ ^11^2 _ 827 , _ 1127 . ^ 91 . ^ 73 .A^ 



04 = — cPA - —a^ + 6a* - 3a^ A + _a^ cf" a* + — a^ + — a'^ A 

4 V 9 / V 3 27 36 4 16 / 

7 5-5 /497_4 ,^ , ^ f:\^ 21 g--6 

a5 = --a''A - ( — o^ + lOa^ -5aMA ag = -"3"^ -^ 
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